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We have studied the wetting properties of *He adsorbed on the surface of heavy alkali metals by 
using a non-local free-energy density-functional which describes accurately the surface properties of 
liquid ''He in the temperature range < T < 3 K . Our results for ''He on the Cs surface give both 
the temperature dependence of the contact angle and the wetting temperature in good agreement 
with the experimental findings. For the *He/Rb system we find that a wetting transition on the Rb 
surface occurs at T ~ lAK, whereas the experiments show either wetting down to T=0 or a wetting 
transition at T ~ 0.3 K. We suggest that this disagreement is due either to an inaccuracy of the 
fluid-substrate potential used in our calculations or the consequence of substrate roughness, which 
is known to affect the Rb surface, and whose effect would be to lower the wetting temperature. The 
sensitivity of the wetting transition to the He-surface potential are stressed for the He/Rb surface, 
which may justify the controversial experimental findings for this system. 



PACS: 68.10.-m , 68.45.-v , 68.45. Gd 
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I. INTRODUCTION 



The existence and nature of the wetting transitionEl at a temperature above the triple point Tt (and below the 
critical temperature) is known to be the result of a delicate balanceJpetween the interatomic potential acting among 
the atoms in the liquid and the atom-substrate interaction potentials. When the latter dominates, the fluid tends to 
wet the surface. Because of the extremely weak interaction between He atoms, liquid "^He in contact with almost any 
substrate spreads to form a continuous film over the surface, so that vapor and substrate are never in contact, thus 
making "^He the ideal wetting agent. 

A remarkable exception to this behavior is found when '*He is adsorbed on heavy alkali metal surfaces: due to the 
large electron spill-out at their surface, alkali, metals provide the weakest adsorption potentials in nature for He atoms. 
Based on this observation, it was suggestedHa that wetting and prewetting transitions of '''He ought to occurr on Cs 
(which exerts the weakest attraction to He atoms) and, perhaps, Rb substrates. Actually, experimental evidence that 
He* on a Cs substrate exhibits a wetting transitionQ^B at « 2 K has been collected. The results for ''He op Rb are 
more controversial, since both a prewetting transition accopapanied by complete wetting down to T = Ko, (this is 
analogous to triple point wetting) and a wetting transitionllJ at a nonzero k, 0.3 K ha,ve been reported. On the 
more attractive substrate of potassium, wetting behavior has been observed for all T > OM. 

Numerical simulations have proven to be a very useful tool to study the wetting properties of fluids. In particular. 
Density Functional (DF) methodst3 have become increasingly popular in recent years because of their ability to 
describe inhomogeneous fluids and phase equilibria. A comparison with "exact" Monte Carlo simulation results shows 
that, once the long-range (van der Waals) attractive forces exerted by a surface on the gas atoms are included in 
the free energy functional. Density Functional methods provide a qualitatively (and most often quantitatively) good 
description of the thermodynamics of gas adsorption on a solid surface and correctly predict a large variety of phase 
transitions (wetting, prewetting, layering, etc.). 

In particular, a quite accurate description of the T = Op*roperties of liquid '^He has been obtained within a DF 
approach by using the energy functional proposed in RefJ13 and later improved in Ref.lU. This phenomenological 
functional has been jiddeljt tested in a variety of problems involving inhomogeneous phases of liquid *He (surfaces, 
droplets, films, etc.)Qoil3, and is known to give a good description of ^^^tfi-fT = properties of the liquid-vapor 
interface and also of some properties of liquid *He on various substratesES-Ej (for a thorough comparison between 
these and other functionals used to describe liquid *He at T = 0, see RefrJ). A very remarkable example of success 
of such theoretical framework is the prediction of the non-wetting behavior of liquid *He on heavy-alkali metals at 
zero absolute temperaturea. r~\r~\ 

In the present work we extend the theory of RefsJl^^o to finite temperatures and apply the functional to the study 
of wettinw)roperties of *He adsorbed on alkali metal surfaces. A preliminary account of this work has been presented 
elsewhereO. 

A prp/ious attempt to construct a semi-empirical finite-temperature Density Functional for "^He has been reported 
in Ref.c3. There, a local functional of the density was used (i.e. where the He-He interaction terms correspond to 
an attractive two-body contact force), with a number of parameters adjusted to guarantee the correct liquid- vapor 
coexistence. The agreement with the experimental surface tension of liquid *He at saturated vapor pressure was also 
imposed by construction. This functiopal has been applied to the study of the behavior of the liquid- vapor interface 
thickness as a function of temperatureE3. 

It is well knownll3, however, that a local theory gives incorrect fluid structures in situations where the liquid is 
perturbed on a microscopic scale, as itiiappens close to a solid substrate. This failure is due to the zero-range nature 
of the Skyrme interaction used in Ref.E3, which completely neglects the short-range correlations associated with the 
interatomic potential hard-core. 

For this reason, in our work we rather use as a starting point for our finite-temperature calculations the non- 
local form of the Density Functional proposed in Refo, where a more xcalistic finite-range interactions between He 
atoms is included. Contrary to what happens with the local functionalEj, this non-local functional with parameters 
fitted to bulk properties provides the correct asymptotic behavior of the surface potential profile without any ad-hoc 
prescription, which is a desirable feature in order to properly describe the liquid-vapor interface. 

II. METHOD 

In the Density- Functional approach the grand-canonical free-energy is considered a functional of the density p{r): 

nHe[p]=F[p]-H j dv p{v) (1) 
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where fi is the chemical potential and F[p\ is the free-energy density functional. 
We propose to write F[p] as: 



fdv (2) 



where 



/ = f quant [p] + fmt [p] (3) 



The first term in the r.h.s. takes into account the quantum kinetic energy and the free energy density of an ideal 
Bose gas Ji^: 

! quant = -^(^Vpf + C4/id (4) 

with C4 an adjustable parameters tcy-he determined by the fitting procedure detailed below. 

As in the original T=0 functionaltil, fint accounts for He-He interactions (in a Hartree-like fashion) and for corre- 
lations effects: 

/»*[p(r)]-^ / dr'p(r)p(r')F,(|r-r'|) 
+ |p(r)(p-r)2-|-|p(r)(p-r)3 

-^as [ dr'Gi\r - r'|)(l - ^)Vp(r)Vp(r')(l - ^). (5) 
4M J pos pos 

The first term contains a two-body He-He pair potential Vi{r) screened at distances shorter than a characteristic 
length hi {Vi — for r < hi), while the second and third terms contain pr, i.e. the average of the density over a 
sphere of radius hi , and account for the increasing contribution of the hard-core He- He repulsion when the density 
is increased. The last term contains the gradient of the density at different points and corresponds to a non-local 
correction to the kinetic energy. 

We remark at this point that the quantum kinetic energy term appearing in (^) (as well as in other phenomenological 
density functional, theories used to describe inhomogeneous "^He) is derived from the total density p{r). However, as 
discussed in RefJHJ, a complete Density Functional theory of Bosc-condcnsed liquids would rather involve functionals 
of both the local density p{r) and the local order parameter $(r) = ^ Pc(r), Pc(r) being the local condensate density. 

In particular, the counterpart of the first term in (^) would be proportional to (Vy^)^. A rationalization of why 
the simplified form of the density functional used here (and in all other applications of DF theory to liquid ''He as 
well) leads to a satisfactory description of the surface can be traced to the fact that the two quantum kinetic eneijgy 
terms coincide in the low-density surface region, which is important for the interface problems, where pdr) ~ p(r)E3. 

In the T=0 functional, the free parameters hi, C2 and C3 were adjusted in order to reproduce the experimental 
values of the density, of the energy per atom, and of the compressibility for bulk liquid ^He at zero pressure, while 
the width Iq of the gaussian function G and the parameters as and pos in Eq. (^) were adjusted to reproduce the 
overall shape of the experimental static,*f!sponse function of bulk liquid '*He. For a detailed description of the various 
terms in (||) we refer the reader to Ref.tJ. 

We follow the same philosophy at finite temperatures, i.e. we consider the free energy density of a uniform system 
(/3=costant) as given in Eqs. (Q,^) : 

f^F/V= \hp^ + ic2p3 + ic3p4 + C4/.,(p, T) (6) 

where b = J dr Vi{r)—4:Tr di" r'^Vi{^), and minimize the grand potential with respect to p to get the chemical 
potential: 

pip, ^) = " + ^^'^P^ + ^"^3^^ + '^'iAiid (7) 



and the pressure 



P{p, T) = pp-f= ^hp" + C2P^ + C3P* + CiP,d (8) 
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If py and pi are the densities of vapor and liquid ^He at saturated vapor pressure at a given temperature T, than 
at coexistence the equiUbrium conditions: 



p{pi,T) = p{pv,T) 



(9) 



and 



P{puT) ^ P{py,T) 



(10) 




(11) 



we get three equations relating the four adjustable coefficients C2, C3, C4, and h. 

As a fourth equation necessary to determine them completely at any temperature of interest we impose that the 
calculated isothermal compressibility, 1/K = p{dP/dp), reproduces its corresponding experimental value: 



The fit has been carried out in the region < T < 3 K, which is the most interesting for the wetting properties of 
"^He we want to investigate. Once b{T) has been obtained by solving the system of four equations described above, 
a straightforward calculations provides hi(T) so that all coefficients entering the functional (H) are determined. The 
calculated coefficients are shown, for some selected temperatures, in Table (|). 

We wish to stress at this point that all unknown coefficients entering the free energy Density Functional of pure 
*He are fitted to bulk properties, i.e. to properties of uniform systems. 

To include the interaction with the substrate, we use a binding potential Vs{z) which describes the interaction 
between the alkali metal, occupying the half space z < 0, and one ^He atom located at a distance z above the ideally 
flat surface. _ . 

The physisorption potential Vs{z) is taken in the form originally proposed in Ref.EEl, i.e. as a sum of a Hartree-Fock 
repulsion and a van der Waals attraction. Both terms are parametrized with a set of coefficients. Very reliable values 
for the parameters entering Vs{z) and describing the interaction of ''He with alkali substrates have been calculated 
from first principle£3. Thege potential have been used to calculate the wetting properties of rare gases (other than 
He) on alkali metal surfacesc2 and also the T=0 profile of a droplet on the Cs surfaceEj, and found to give predictions 
in reasonable agreement with experiments for all these systems. 

The total free energy functional for ^He interacting with the substrate is thus: 



According to Density Functional theory, the equilibrium density profile p{r) of the fluid in the presence of the 
substrate can be determined by applying the variational principle : 



K{pi,T)^K--P'{T) 



(12) 




(13) 



Sn/6p{r) = 



(14) 



to get the Euler-Lagrangian equation: 



{-4?^' + u[pir)] + v;(z)} VpW = mVpW, 



(15) 



where the effective potential U is deflned as C/[/3(r)] = 5F[p{r)]/S^/p and the value of the chemical potential ^ is fixed 
by the normalization condition 




(16) 



N being the total number of ''He atoms. 



4 



III. RESULTS AND DISCUSSION 



Wetting or non wetting is determined by surface tension balancing. The intersection of a macroscopic liquid droplet 
with a solid planar substrate can be characterized by the contact angle @: it goes to zero as T approaches T^, and 
remains zero at all higher temperatures, where the liquid wets the substrate. Below the wetting temperature, the 
contact angle is determined by balancing the forces acting along the contact line and depends on the interfacial 
tensions aij between each pair of coexisting phases through Young-Dupre's equation 

coskd ~ . (17) 

The subscripts I, v and s identify the liquid, vapor and solid, respectively. 

We have calculated the interfacial tensions CTy entering Eq. at selected values of T by using the density 

functional method described in the previous section. 

We first studied the liquid-vapor interface of ^He by solving Eq. ( [l5| ) (with Vg = 0) with the appropriate boundary 
conditions, corresponding to a planar liquid film in equilibrium with its own vapor. In this case the solution depends 
only on the coordinate z normal to the surface, and provides the equilibrium density profile p{z). Our calculations are 
performed within a region of length Zm- Rather large values of Zm and of the film thickness tf are necessary for well 
converged calculations (typically Zm ~ 200 A and tf ^ 150 A), in order to let the system spontaneously reach (i) the 
liquid bulk density in the interior of the film and (ii) the vapor bulk density in the region outside the film. In Fig.l 
we show some selected density profiles close to the liquid-vapor interface. A detailed analysis of the resulting density 
profiles shows that the widiJi-of the interface increases monotonically with temperature. A comparison of our results 
with those obtained in Ref.t3 indicates also that the thickness of the liquicLaiapor interface, at a given temperature, 
is sligthly smaller than that obtained from the zero-range functional of Ref.t3. 

From the calculated equilibrium density profiles one can directly compute the liquid- vapor surface tension (T/„ , from 
the definition a = {fl + PV)/A. Here P is the saturated vapor pressure at temperature T, V is the volume of the 
system and A is the surface area. For the one-dimensional problem considered here, one can write: 

a = (l/2)[ / " f{p{z))dz - M / " p{z)dz + Pz^] (18) 
Jo Jq 

The pressure P can be conveniently calculated from the definition P = (/x/9„ — f{pv)) = {ppi ~ f{pi)), Pv and pi 
being the experimental densities for bulk liquid and vapor, respectively. The chemical potential p is an output of our 
calculation, being fixed by the areal density rij, = N/A = /q™ p{z)dz. Finally, a factor 1/2 appears in the previous 
equation to account for the two free surfaces delimiting the liquid film in our "slab" calculations. 

The dependence of the liquid-vapor surface tension of "^He on temperature has been extensively investigated in the 
past. Surprisingly, the absolute valuej-of pOj^ at T = remains poorly known: predictions from different groups give 
values for aiv{0) differing by up to 69ccHi^. Our calculated values are shown in Fig. 2, where they are compared with 
the available experimental data. The overall agreement is quite satisfactory, given the fact that our functional has 
been fitted to reproduce ^He bulk properties only. Note also that the kink (barely) visible in the experimental data 
at the temperature corresponding to the A-point T\ = 2A7 K appears also in our calculated aiy and reflects similar 
kinks appearing close to T\ in the bulk quantities (density, chemical potential and compressibility) used to fit the 
coefficients entering the free energy (||). 

We next turn to the problem of liquid ^He in the presence of a solid surface. The two additional surface tensions 
asv and agi required to obtain the contact angle O from Eq. (^^ are calculated within our density functional 
approximation by using different boundary conditions from those used to calculate aiy . We confine the fluid between 
two parallel surfaces, separated by a distance large enough to avoid any compression effect. The two surfaces act 
on the fluid with the external potential Vs{z) appropriate to the alkali metal substrate. For selected values of the 
temperature T, we get the structure of the solid-vapor and solid-liquid interfaces by minimizing the grand-potential 
( p^ subject to the additional constraint that the density in a large region between the two surfaces (where bulk 
behavior is expected) is equal respectively to py and pi. The corresponding surface excess energies per unit area 
a = + PV)/{2A) are the solid-vapor asv and the solid-liquid asi surface tensions. Of course the distance between 
the two confining surfaces must be such that the calculated values of asv and asi do not change by further increasing 
their separation. We show in Fig. 3 the equilibrium density profiles of ^He confined between two Rb surfaces, at 
T = lAK. The solid line shows the density for the liquid phase confined between the surfaces, the dotted line (which 
is rescaled for clarity by a factor 30) shows instead the vapor density profile. Note the oscillations of the liquid density 
close to the surface, due to close-packing effects. We have checked, by means of similar calculations with different cell 
sizes, that at least 150 A must separate the two attractive surfaces so that the two fluid-solid interfaces are effectively 
decoupled and the calculation gives a converged value for the surface tensions asi and asv ■ 



5 



According to Eq. (0), wetting of a surface occurs at a temperature such that 

Csv ~ CTsl — CTly (19) 

In Fig. 4 we show the behavior of the two sides of Eq.(^9|) as a function of temperature, for three different alkaH 
surfaces. A common feature of thCpGurves for asv — <^si shown in Fig. 4 is their rather weak temperature dependence, 
as indeed observed in experiments^ The predicted wetting temperature for Cs, ~ 2.1 K, is in good agreement 
with the experirriental value ^ 2 K. Also in the case of '*He on K our calculated results shown in Fig. 4 agree 
with experimentga, indicating wetting at any T > 0. On the contrary, for Rb our calculations predict a wetting 
temperature ^ 1.4 K, in marked disagreement with the expjerimental findings which show instead either complete 
wetting down to T = 00, or a wetting transition at « 0.3 KllJ. 

We show in Fig. 5 the temperature dependence of the contact angle for the ''He/Cs system as obtained from 
our calculations, together with the available experimental data. Although the equilibrium contact angle is a well 
defined thermodynamic quantity, measurements of this quantity invariably show a strong hysteretic behavior, i.e. 
the measured value of the contact angle depends on whether the contact line is advancing or receding. For relatively 
homogenous surfaces, however, it has been shown that the advancing contact angle is a good measure of the equilibrium 
contact angle, so we compare our results with the experimental values for that quantity. The cogtact-angle for. "^He . 
on Cs surfaces has been measured by different groups in recent years by using different techniquesEiEJEa. In RefE^a 
the cornet angle is measured by direct visual inspection of macroscopic droplets of superfluid ^He on a Cs surfaee. 
In Ref.E^ the contact angle is obtained by means of direct optical measurement using interferometric techniquesES, 
whereas in Ref.u it is measured, much more indirectly, from the reduction in pressure due to capillarity, when an array 
of parallel tungsten plates coated with Cesium are immersed into liquid "^He. Our resultSpare-compared in Fig. 5 with 
the aforementioned experimental data. The overall agreement with the results from Ref.E^TtJ is satisfying, although 
a strong discrepancy exist with the data of Ref.El. 

At variance with the overall satisfactory quantitative agreement we get for the ^He/Cs system, our predictions 
seems to be wrong for the ''He/Rb system. We believe that this disagreement is probably due to the fact that our 
assumptions of a rigid and planar substrate are far from being satisfied in this case. Impurities or roughness, in fact, 
are always present on the surface, as indicated by the strong hysterectic behaviour of the contact angle observed in 
all experiments when studying advancing or receding contact lines. This condition may be particularly severe in the 
case of the Rb surface, where defects and roughness are known to be present to a larger extent than on the surface of 
Cs. 

The heterogeneity of the substrate is known to promote wetting, i.e. leads to values of a^y — a^i higher than for an 
ideal substrate. In the case of ^He on Rb, as apparent from Fig. 4, an increase of asv — o'si would shift the wetting 
temperature towards lower values, as indeed found in the experiments. 

The characterization of the wetting behavior of superfluid Helium on rough surfaces is a very difficult task. The 
quality of the substrates used in the experiments cannot be determined at lengthscales smaller than optical wave- 
lengths. The theoretical description of the effecte,pf heterogeneities is also challenging. Only few theoretical results 
are available on the effect of microscopic disordeii230 on the wetting properties of fluids adsorbed on weakly attractive 
surfaces. 

The effect of model disorder (on an atomic scale) on the wetting properties of Ne on the Mg surface has been studied 
by Curtarolo et alE3. One interesting conclusion contained in this paper is that roughness destroy the discontinuous 
nature ofJJae wetting transition, even if a trace of discontinuities remain in the calculated isotherms. 

In Ref.LJ a Lennard-Jones fluid and its wetting properties on molecularly rough surfaces have been studied by 
using Molecular Dynamics. Contrary to the conventional belief that surface roughness reduces the contact angle thus 
making it easier to wet a rough surface than a smooth one, the surprising conclusion of this work is that the contact 
angle is larger for the rough surface than for the smooth one. 

In order to have further insight into this important issue, we are currently investigating the effect of a realistic 
distribution of microscopif; disorder on the wetting transition of ^He on the Rb surface within the Density Functional 
approach described hcrtH. 

An alternative explanation for the disagreement with the experiments for the ^He/Rb system is that the He-Rb 
surface potential used in our calculations is inaccurate. Our overestimate of the wetting temperature would be a sign, 
in this case, that the calculated potential is slightly less attractive than required to provide a result closer to the 
experimental one. 

We have thus tried a modified He/Rb surface interaction by changing sligtly the two parameters entering the 
repulsive part of the ab-initio potential Vs{z) describing the He-surface interaction. The modified potential, has a 
minimum of depth D — 0.661 meV located at Zmin = 8.46 oq, to be compared with the values D — 0.629 meV 
and Zmin — 8.60ao obtained with the original parameters. Such a tiny variation in the adsorption potential (^ 5 % 
increase in the potential well depth) has quite a large effect on the He/Rb wetting temperature, as can be judged from 
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the curve labeled Rb* in Fig. 4. With this modified potential the wetting transition moves very close to T = 0, as in 
the experiments. Interestingly, in this case the curve —asi(T) almost coincides with the curve aiv{T) for a range of 
T values, between and 0.3 — 0.4 K (at such temperatures the contribution agv is negligible in our approximations). 
This could provide a reason for the different experimental estimates of the wetting temperature (either T = or 
T ~ 0.3 — 0.4 K) by different groups: in the region where the two curves almost overlap their slope is very small, and 
thus small variations in the quality of the substrate may result in rather different wetting temperatures. Moreover, 
this could also suggest a different character of the wetting transition for He/Rb, since similar slopes in —asi{T) and 
ai-uiT) imply a small discontinuity in dcos{Q)/dT, i.e. a quasi-continuous transition. 

IV. SUMMARY 

We have extended to finite-temperatures the non-local Density Functional proposed in recent years to describe 
the *He properties at T^OGja, and studied the wetting behavior of ^Hc adsorbed on alkali metal surfaces. The 
Density Functional depends on a number of temperature-dependent phenomenological parameters which are adjusted 
to reproduce bulk experimental properties of liquid ^He at saturated vapor pressure. We find that the resulting 
functional describes accurately the properties of the liquid-vapor interface of the free ^He surface. In the presence 
of an external potential simulating a planar, solid surface whose binding properties reproduce those of alkali metal 
surfaces, it also accurately reproduces the properties of the ^He/Cs system, giving the temperature dependence of the 
contact angle in good agreement with experiments. The contact angle vanishes at T ~ 2.1 K, in agreement with the 
experimental measure for the wetting temperature, T ^ 2K. For ^He on a K surface, we find wetting at all nonzero 
temperatures, as found also in experiments. For ^He on the Rb surface, however, our calculations fail to reproduce 
quantitatively the experiments, indicating a wetting transition at T ^ lAK . The discrepancy with experiments could 
be either due to our neglect of surface inhomogeneities, which are known to be present on this surface to a larger 
extent than on a Cs surface, or to some inaccuracy in the theoretical determination of the surface-adatom potential 
used in our calculations. If this is the case, our results indicate that the correct potential should be slightly more 
attractive than currently believed. 
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Figure captions: 



Figure 1: Calculated "^He density profile near the liquid- vapor interface. Solid line: T ^ OK; dotted line: T — 1.5 K; 
Dashed line: T ^2K; Dash-dotted line: T = 2.5 K. 

Figure 2: Surface tension ofJicwid, *Hej— Dots: calculated cr;„; Squares, diamonds, triangles and crosses show the 
experimental results from Ref.E3,E3,Eil andc2l, respectively. Lines are only a guide to the eye. 

Figure 3: Equilibrium density profiles used to calculate a^i (upper curve, with solid line), and asv (lower curve, 
with dotted line), respectively. 

Figure 4: Calculated surface tensions for the ^He/alkali systems. Squares: ai^; Dots, diamonds and crosses show 
{asv — cfsi) for '*He/Cs, *He/Rb and ''He/K, respectively. The set of points labeled Rb* are obtained by using a 
modified ^He/Rb surface interaction, as explained in the text. Lines are only a guide to the eye. 

Figure 5: Temperature dependence otitlw; ^Hoj/Cs contact angle: squares: this work; triangles, dots and crosses 




8 



^ J.W.Cahn, J.Chem.Phys. 66, 3667 (1977). 

^ See, for example, S. Dietrich in Phase transitions and critical phenomena vol. 12, cd. by C. Domb and J. Lebowitz, Academic 

Press, San Diego (1988); B. Evans and M.H.W. Chan, Phys. World, 9 , 48 (1996). 
3 E.Cheng, M.W.Cole, W.F.Saam and J.Treiner, Phys. Rev. Lett. 67, 1007 (1991). 
* E.Cheng, M. W. Cole, W. F. Saam, and J.Treiner, Phys. Rev. 46, 13967 (1992). 

^5 P. Taborek and J.E. Rutledge, Phys. Rev. Lett. 68, 2184 (1992); J.E.Rutledge and P.Taborek, Phys. Rev. Lett 69, 937 

(1992). 

® P.J.Nachcr and J.Dupont-Roc, Phys. Rev.Lett 67,2966 (1991). 

^ K.S. Ketola, S. Wang, and R.B. Hallock, Phys. Rev. Lett. 68 201 (1992); R. B. Hallock, J. Low temp. Phys. 101, 31 (1995). 
® J. Klier, P. Stefanyi and A.F.G. Wyatt, Phys. Rev. Lett. 75, 3709 (1995);J.Klier and A.F.G.Wyatt, J.Low.Temp. Phys. 110, 

919 (1998). 

^ J. A. Phillips, D. Ross, P. Taborek, and J. E. Rutledge, Phys. Rev. B 58, 3361 (1998). 

1° A. F. G. Wyatt, J. Klier and P. Stefanyi, Phys. Rev. Lett. 74 1151 (1994); J. Klier and A.F.G Wyatt J. Low Temp. Phys. 
113, 817 (1998). 

" J. A. Phillips, P. Taborek, and J. E. Rutledge, J. Low Temp. Phys. 113, 829 (1998). 
" R.Evans, Adv. Phys. 28, 144 (1979). 

" J.Dupont-Roc, M.Himbert, N.Pavloff and J.Treiner, J. Low Temp. Phys. 81, 31 (1990). 

" F.Dalfovo, A.Lastri, L.Pricaupenko, S.Stringari and J.Treiner, Phys. Rev. B 52, 1193 (1995). 

^^F.Dalfovo, Z. Phys. D 29, 61 (1994); M.Casas, F.Dalfovo, A.Lastri, Ll.Serra, and S.Stringari, Z. Phys. D 35, 
67 (1995); A.Lastri, F.Dalfovo, L.Pitaevskii and S.Stringari, J. Low Temp. Phys. 98, 227 (1995); F.Dalfovo, 
A.Franchctti,A.Lastri,L.Pitacvskii and S.Stringari, Phys. Rev. Lett. 75, 2510 (1995);J. Low Temp. Phys. 104, 367 (1996). 
M.Barranco and E.S.Hernandez, Phys. Rev. B 49, 12078 (1994). 

" F.Ancilotto, E.Cheng, M.W.Cole and F.Toigo, Z. Phys. B 98, 323 (1995). 
F.Ancilotto and F.Toigo, Phys. Rev. B 50, 12820, (1994). 

" F. Ancilotto, A.M. Sartori and F. Toigo, Phys. Rev. B 58, 5085 (1998). 

2° L.Szybisz, Eur. Phys. J. B 14, 733 (2000). 

F. Ancilotto, F. Faccin and F. Toigo, Proceedings of the 14th Symposium of Thermophysical Properties (Boulder, Colorado, 
June 25-30, 2000), to be published in Fluid Phase Equilibria. 

A.Guirao, M.CentcUcs, M.Barranco, M.Pi, A.Polls and X.Vinas, J. Phys. Cond. Matter. 4, 667 (1992). 

R.K.Pathria, Statistical Mechanics (Oxford, Pergamon, 1972). 
2^ A. Griffin and S.Stringari, Phys. Rev. Lett. 76, 259 (1996). 
'^^ J.Maynard, Phys. Rev. B 14, 3868 (1976). 
2® E.Zaremba and W.Kohn, Phys. Rev. B 15, 1769 (1977). 

^'^ A. Chizmoshya, M.W.Colc and E. Zaromba, J. Low Temp. Phys. 110, 677 (1998). 
2* F. Ancilotto and F. Toigo , Phys. Rev. B 60, 9019 (1999). 

M.Ino, M.Suzuki and A.J.Ikushima, J. Low Temp. Phys. 61, 155 (1985). 

K.R.Atkins and Y.Narahara, Phys. Rev. 138, A437 (1965). 
3^ K.N.Zinoveva and S.T.Bolarev, Sov. Phys. JETP 29, 585 (1969). 

^"^ P.Rocho, G.Dcvillc,N.J.Appleyard and F.I.B.Williams, J. Low Temp. Phys. 106, 565 (1977). 

E.Rolley and C.Guthmann, J.Low.Temp. Phys. 108, 1 (1997). 

D.Ross, J.E.Rutledge, P.Taborek, Science, 278, 664 (1997). 

J.E.Rutledge, D.Ross, P.Taborek, J. Low Temp. Phys. 113, 811 (1998). 
3« S.Curtarolo, G.Stan, M.W.Cole, M.J.Bojan, and W.A.Steele, Phys. Rev. E 59, 4402 (1999). 

J.Z.Tang and J.G.Harris, J. Chcm. Phys. 103, 8201 (1995). 

S.M. Gatica, F. Ancilotto and F. Toigo, unpubUshed results. 



9 



TABLE I. 



T{K) 


PI (A"') 


Pv (A-') 


b (i^l^) 


C2 (JcA"^) 


C3 (KA^) 


C4 


hi (A) 


0.0 


0.0218354 


0.0 


-719.2435 


-24258.88 


1865257 


0.0 


2.19035 


0.4 


0.0218351 


0.441987(-11) 


-714.2174 


-24566.29 


1873203 


0.98004 


2.18982 


0.6 


0.0218346 


0.452619(-8) 


-705.1319 


-25124.17 


1887707 


0.99915 


2.18887 


0.8 


0.0218331 


0.133640(-6) 


-690.4745 


-26027.12 


1911283 


0.99548 


2.18735 


1.2 


0.0218298 


0.495337(-5) 


-646.5135 


-28582.81 


1973737 


0.99666 


2.18287 


1.4 


0.0218332 


0.147990(-4) 


-625.8123 


-29434.03 


1984068 


0.99829 


2.18080 


1.6 


0.0218453 


0.346234(-4) 


-605.9788 


-30025.96 


1980898 


1.00087 


2.17885 


1.8 


0.0218703 


0.684096(-4) 


-593.8289 


-29807.56 


1945685 


1.00443 


2.17766 


2.0 


0.0219153 


0.119400(-3) 


-600.8313 


-27850.96 


1847407 


1.00919 


2.17834 


2.1 


0.0219500 


0.151732(-3) 


-620.9129 


-25418.15 


1747494 


1.01156 


2.18032 


2.2 


0.0219859 


0.188202(-3) 


-619.2016 


-25096.68 


1720802 


1.01436 


2.18015 


2.4 


0.0218748 


0.275044(-3) 


-609.0757 


-26009.98 


1747943 


1.02130 


2.17915 


2.6 


0.0217135 


0.383550(-3) 


-634.0664 


-23790.66 


1670707 


1.02770 


2.18162 


2.8 


0.0215090 


0.516536(-3) 


-663.9942 


-21046.37 


1574611 


1.03429 


2.18463 


3.0 


0.0212593 


0.676814(-3) 


-673.6543 


-20022.76 


1535887 


1.04271 


2.18562 
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